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Navigating networks comprising the interdependencies between the elements of large complex systems and 
databases has been shown as a useful tool for ranking the importance that these constituents have on their global 
organization. A key example is Pagerank navigation which is at the core of the most used search engine of the 
World Wide Web. With the advent of quantum technologies we have learnt that some classical tasks can be out- 
performed when quantum rules enter into play. Here we show that it is possible to define a quantum navigation 
providing a unique ranking reflecting the network topology. We demonstrate that the computation of this quan- 
tum ranking is faster and solves degeneracies of classical ranks. Moreover, we show that quantum coherence 
captures a more global view of the interdependencies between elements, highlighting new hierarchical features. 
Contrary to the usual quantum algorithms, our quantum rank does not need a quantum computer, and we show 
its implementation in real networks. 
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I. INTRODUCTION 

The search for information in the World Wide Web 
(WWW) through search engines has turned into a daily habit 
and an essential tool to fulfill most of our work duties. An 
ideal search engine looks for the information the user is query- 
ing amongst billions of webpages in real time, and produces 
a ranking of the results sorted according the user expecta- 
tions. Although not being among the first search engines 
available, the Google search engine was the first to achieve 
these goals efficiently, establishing one of the milestones of 
the digital era. Its main novelty was to classify and rank 
webpages based on the interrelations created between them 
through the hyperlinks iH, rather than using only their intrin- 
sic features (such as the page content). Google's ranking al- 
gorithm, known as PageranklQ] (PR), is rooted in a diffusion 
process that mimics the user's navigation through webpages 
as the motion of a random walker following hyperlink path- 
ways. 

After Google's global success, multiple applications of PR 
navigation have blossomed whenever a large dataset can be 
mapped into a complex network encoding the interdependen- 
cies between items. Some examples include the classification 
of species within ecosystemsist] or the evaluation of scientists' 
impact in different research disciplines |4{|. In these examples 
PR has successfully produced a classification in which the im- 
portance of each element accounts for its status within the 
whole system. Recent studies[5] have shown that what un- 
derpins the success of PR navigation in the classification of 
network elements according to theirglobal role is the scale- 
free nature of most real networks ja-fS!. highlighting the influ- 
ence that the structural properties of real networks have on the 
collective outcome of the dynamical processes taking place on 
top of themigjjisll. 

Very recently, the dynamical setting of network-based dif- 
fusion has been extended to the quantum domain fl6ll . In this 
framework, quantum random walks have already shown their 



potential for practical purposes as they offer a quadratic speed 
up with respect to the classical algorithms for searching in un- 
sorted databasis lflTlflsll . Considering these recent advances, 
it is tempting to explore the potential application of quan- 
tum random walks for ranking elements in large complex net- 
works, i.e., a quantum rank (QR) algorithm, that improves PR 
based on the benefits that quatumness may introduce lfl9ll . To 
achieve this challenge we need to solve two important issues, 

( i) how to extend the quantum formulation of a walk in a sim- 
ilar way as PR did with classical walks so to find a reliable 
and unique ranking reflecting the topology of the network and 

( ii) how to identify the role that quantum coherence of the QR 
dynamics plays to capture new features about the organization 
of large complex networks. 



In the present work we show how to design a QR navi- 
gation as an hybrid classical-quantum walk. Our theoretical 
approach relies on a master equation for the walkers motion 
interpolating between pur ely quantum coherent dynamics and 
the classical diffusion ||20|] . We explicitly show that a suitable 
choice of the master equation for the walkers motion allows 
to find a stationary solution yielding a unique and reliable 
QR. We show that the interplay of quantum coherence with 
the classical hopping dynamics yields an optimal operation 
point so that QR navigation becomes remarkably faster than 
its classical counterpart. We will further analyze the effects 
that quantum coherence produces by comparing the rankings 
produced by QR and PR navigations. As a result we show 
that QR is able to split the degeneracies found in the ranking 
of PR and, more importantly, to unveil hidden hierarchies in 
the rank of lowly connected items that scape the resolution of 
PR. 
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II. RESULTS 

A. Quantum versus Classical Random Walks in complex 
networks. 

The architecture of a complex network is described by 
means of a collection of N nodes and L edges accounting for 
the elements of the system and their pairwise interactions re- 
spectively. These latter interactions are encoded into a N x N 
adjacency matrix A = {aij}, so that Uij = 1 when a link 
from i to j exists while = otherwise. Given this simple 
formulation one can design a variety of ways for exploring the 
network topology, being a random walk the most simple way 
to navigate networks. 

A random walk is usually described as a time-discrete pro- 
cess. At each time step, a walker hops from a node i to another 
node j provided such connection exists, i.e., aij = 1. The 
dynamics can be formulated in terms of a transition matrix 
n — {TTij} being the entry TTy the probability that the walker 
jumps from itoj. For usual random walks the transition ma- 
trix n is defined as 11^^ = aij/ki, where ki ~ X^jLi '^ij is 
the out-degree of the node i, i.e., the number of links outgoing 
from node i. The dynamics of the random walk is monitored 
by the time evolution of a vector = {pi} whose i-th com- 
ponent accounts for the probability that the walker is placed 
at node i at time t. Then, starting from some initial occupa- 
tion probability, p", the dynamics of the walker follows the 
time-discrete rate equation; 

JV 

^'* = E^y■P^'' withi = i,...,iv. (1) 

The time-continuum version of a classical random walk can 
be also easily obtained from the above rate equation as: 

dtp = {ii-i)p, (2) 

where I is the identity matrix. On both cases, the itera- 
tion of Eq. ([T]i and Eq. ^ ends up in a stationary state 
p* = limt^ooP*- This stationary solution p* is the cor- 
nerstone of any diffusion-based ranking as each of its com- 
ponents p* accounts for the importance of the corresponding 
node. 

The usual map between classical and quantum random 
walks consists in identifying the jump probabilities (contained 
in n in the classical formulation) with ladder operators of a 
tight-binding Hamiltonian whereas walkers are described by 
two level systems, i.e., qubits lfl6ll . Quantum walks have been 
discussed in both continuum and discrete versions ones. How- 
ever, since they were originally conceived as implementations 
of quantum algorithms, all of them share an important prob- 
lem regarding the design of a QR: as quantum walks are re- 
versible, they do not have an stationary solution for the occu- 
pation probabilities, thus making the definition of QR in these 
schemes quite challenging 12 10 (see Appendix for a deeper dis- 
cussion). 

Here we take a different route. In a closer analogy to the 
classical Markovian dynamics for the occupation probabilities 



defined in Eq. (O, we write its quantum analogue consider- 
ing the N X N quantum density matrix g. In this case, the 
map between quantum and classical random walks lies on the 
fact that the diagonal elements of g account for the occupation 
probabilities of each node; = gu = Pi- Following this 

prescription, it turns out that any Markovian quantum evolu- 
tion can be written in the form ll22ll (see Appendix); 

^=-i(l-a)[i/,e] (3) 

where H is the Hamiltonian {H = W) incorporating the 
(quantum) coherent dynamics, whereas 7(ij) > and the 
operators j)} (the so-called dissipators) are the ones re- 
sponsible of irreversibility. The parameter a £ [0,1] quanti- 
fies the interplay between unitary and irreversible dynamics. 
The above formulation avoids the aforementioned problem as- 
sociated to unitary evolution, i.e., the absence of a stationary 
distribution for the quantum occupation probabilities. More- 
over, the irreversible part of Eq. (O can accommodate any 
classical random walk by choosing the definition of the dissi- 
pators. For instance, the usual random walk can be incorpo- 
rated by choosing; = and7(j j) = n^j. In this 
case, in the limit a ~ 1, the equations for the diagonal ele- 
ments, gjj = pj meet those of the usual classical walk, Eq. 
(|2]i. Thus, Eq. ([3]) interpolates between classical (a = 1) and 
purely quantum random walks (a = 0). Between both lim- 
its the dynamics incorporates the interplay between classical 
hopping and quantum coherence thus making navigation sen- 
sitive to non-local correlations between the nodes and paving 
the way to a more global ranking. 

B. Quantum rank. 

As mentioned previously, once the network navigation is 
defined the ranking of its nodes is done by evaluating each of 
the components of the stationary distribution, p*. This sta- 
tionary distribution corresponds to the eigenvector of the tran- 
sition matrix 11 with (the largest) eigenvalue 1, Hp* = p*. 
However, the existence and unicity of p* is only guaranteed 
provided the transition matrix fulfills the Perron-Frobenius 
theorem, i.e., it is irreducible and aperiodic ll23l . In order to 
satisfy this latter condition in any kind of network topology, 
the PR navigation incorporates a long-distance hopping prob- 
ability to the usual classical walk. This ingredient transforms 
the usual transition matrix, 11, into the so-called Google ma- 
trix; 

G = qU+{l-q)F , (4) 

where q G [0, 1] is a parameter whose value is typically set 
to q = 0.9 and F is the long-distance hopping matrix, Fij = 
1/{N — I)ifiy^j and otherwise. By doing so, the transition 
matrix G has a unique eigenvector with maximal eigenvalue 
1. Thus, the ranking of websites is related to the components 
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of the eigenvector corresponding to the (unique) eigenvalue 1 
of the Google matrix, Eq. (HJi. 

In the quantum case we face a similar problem. The defini- 
tion of the dissipators in Eq. ^ with the form of the classical 
operators, = N)OI and 7(ij) = Hy, does not guaran- 

tee a unique stationary distribution for any network topology. 
However, it is possible to prove (see Appendix) that, when dis- 
sipators are defined as L/^^ j-j = with 7(ij) — Gij, Eq. 
(O has always a unique stationary solution. In a nutshell, the 
patching mechanism used in the Google matrix is generalized 
here to the quantum case in a natural way. As a consequence, 
by incorporating the PR navigation in the irreversible part of 
Eq. (O we can compute a reliable QR that incorporates the 
effects of coherences in the navigation. 



C. Quantum speed up. 

We now study the performance of our QR compared with 
that of PR. First we focus on one important aspect for the com- 
putation of a rank in large networks: the convergence time, r, 
to the stationary state of the navigation dynamics, g*, i.e., the 
rank. Being linear, the time evolution for the master equation 
(O can be written in a compact form as g ^ J0.[g]. Conse- 
quently, the evolution is determined by the spectrum of C and 
the convergence to the stationary state is upper-bounded by 
T = |Rc[Ai] 1^^, with Ai being the first non zero eigenvalue 
ofC. 

To study the convergence time of QR we have considered 
two of the most paradigmatic network topologies: Erdos- 
Renyi (ER) and scale-free (SF) networks. These two networks 
differ in the distribution P{k) of the number of contacts (de- 
gree), k, per node. In ER graphs most of the nodes share 
the same degree and P{k) is shown to be a Poisson distribu- 
tion (see Fig. l.a). However, real-world networks follow a 
SF pattern characterized by a power-law degree distribution: 
P{k) ^ k^^ (see Fig. l.b) where the exponent lies in the in- 
terval 2 < 7 < 3 so that the variability of the degrees found 
in the network becomes extremely large. 

In Figs. [T]I.b and l.d we plot r as a function of the 
classical-quantum interplay included in the QR navigation 
[quantified by a in Eq. (|3]l] for SF and ER networks re- 
spectively. We observe that for a — !> (approaching the 
fully-quantum limit) the value of t diverges, as expected, 
since in this limit the dynamics becomes fully-unitary. On 
the other hand, in the classical limit, a — > 1, the conver- 
gence time tends to its classical counterpart, reaching 
T^^ at a = 1. The striking result appears for intermedi- 
ate values of a: the curve r(a) shows a minimum for some 
< Q!°P* < 1. Thus, there exists an optimal value a"?'* point- 
ing out that an hybrid classical-quantum navigation outper- 
forms significantly the time needed for ranking with PR. This 
enhancement resembles the quantum stochastic resonance- 
like phenomenon ll24ll previously found in other contexts, such 
as entanglement generation 1251 or efficient light-harvesting 
transport 11261? , |2711. In our case, it is the competition be- 
tween the coherent and the dissipator parts what accelerates 
the convergence to the stationary distribution of walkers, i.e., 
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FIG. 1. Convergence times: QR versus PR. In this figure we show 
the evolution of the convergence time to the stationary solution of 
QR, r'^^, as a function of the quantum-classical interplay, a, intro- 
duced in the QR navigation. In the top part of the figure we focus on 
ER graphs in which the probability, P{k), that a node is connected 
to k other nodes (the degree distribution) follows a Poisson distri- 
bution as show in panel (a). Panel (b) shows the evolution r'^^(a) 
compared to the convergence time of PR, r^^, in an ER graph of 
— 200 nodes. It is clear that QR reach faster the stationary so- 
lution (the final rank) for a broad range of a values, reaching the 
maximum outperformance for a"''* ~ 0.54. The architecture of SF 
networks is rather different as shown in their power-law degree dis- 
tribution, see panel (c). In panel (d) we plot the evolution T'^^{a) 
in a SF network with P{k) ~ k~'^ and iV = 200 nodes. We ob- 
serve again that for a range of a values QR converges faster than 
PR to a stationary ranking, reaching the maximum outperformance 
at a"*'' ~ 0.68. The speed up in the converge of QR with respect 
to that of PR is thus observed in both ER and SF networks. In panel 
(e) we show, for SF networks, the scaling of the ratio t-Q^/t^^ at 
a°P* with the size of the network, A'^. As observed, the former ratio 
remains roughly constant, thus showing the robustness of the speed 
up provided by QR. 



the ranking. Finally, in Fig. [T]l.e we show the scaling of the 
ratio I as a function of the size N of the network for 
SF networks. It is clear that the enhancement provided by the 
QR navigation is not a matter of the finiteness of the network 
as it remain roughly unaltered for larger values of A^. 



D. Solving PR degeneracies 

After providing a reliable definition of QR that provides a 
unique classification of nodes, the speed up shown in the pre- 
vious paragraph supports its feasibility as a navigation-based 
ranking method. At this point the important issue relies on 
the novel features introduced in the ranking itself by the in- 
terplay of quantum and classical effects that QR incorporates. 
To perform a comparison between QR and PR we first rely 
on a simple network[|29;l in which a qualitative rank can be 
made at first sight to check the output of different rank algo- 
rithms. In particular, we will test the rank provided by three 
navigation dynamics: (i) the usual random walk, (ii) PR, and 
( Hi) QR. The simple network is a directed graph composed of 
eight nodes (see Fig. 2. a) that can be described as the sum of 
a densely connected core (composed of nodes 1, 2, 3 and 4) 
plus a peripheral cyclic circuit (composed of nodes 5, 6 and 
7). Given the simple structure of this graph, it is easy to have 
some insight about the expected ranking. The nodes belong- 
ing to the densely connected core will accumulate most of the 
walkers' flow in the graph while, in the peripheral part, node 
5 will receive all the outgoing flow from the connected core. 
Therefore, these nodes should occupy the top positions of a 
proper ranking algorithm. 

In Fig. |2f .b we show a table containing the ranks produced 
by the aforementioned three navigation schemes. The three 
ranks show that the four nodes belonging to the densely con- 
nected core obtain the largest ranks, being node 2 the most 
central node as it is the unique element receiving extra-flow 
from the peripheral part of the graph. On the other hand, 
nodes 1 and 4 obtain a similar centrality score highlighting 
their equal role within the graph. The differences between 
the three navigation schemes show up in the classification of 
the peripheral nodes. First, the classical random walk assigns 
similar scores to the couples of nodes (5, 7) and (6, 8). This 
degeneracy is a product of some local similarities between the 
two nodes of each couple. Nodes 8 and 6 allocate the same 
incoming flow of walkers (that coming from node 7) while in 
the case of nodes 7 and 5 it is clear that all the flow reach- 
ing node 5 from the connected core passes directly to node 7. 
However, at a global scale, the role of these four peripheral 
nodes is not identical, thus challenging the reliability of the 
classical walk as a method to capture the differences between 
the global role of nodes. 

The degenerate rank produced by the classical random walk 
demands the use of navigation algorithms with a broader in- 
formation horizon so to unveil the differences in the role 
played by peripheral nodes at a global scale. PR incorporates 
this nonlocal ingredient by means of the long-range hops in- 
corporated in the Google matrix (|4]i. From Fig|2^.b we ob- 
serve the former degeneracies seems to be partially solved by 
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FIG. 2. QR versus classical ranks in a small graph. In panel (a) 
we show a simple directed graph composed of 8 nodes. This graph 
can be briefly described as the sum of a densely connected (com- 
plete) subgraph (composed of nodes 1, 2, 3 and 4) plus a 3-cycIe 
(composed by the closed circuit 5 — >■ 7 — >■ 6 — > 5). These two 
structures are joined together by means of a link 3 — >■ 5 and a di- 
rected path 7 — >■ 8 — !• 2. In panel (b) we show the rankings obtained 
via two classical navigation schemes, the usual random walk (RW) 
and PR with q = 0.9, and our QR. For the three navigation-based 
rankings the first 4 nodes (composing the densely connected subset) 
are the most central ones, being headed by node 2 (as it is the only 
one receiving extra flux from the peripheral part) and followed by 
node 3 (being the only one with one outgoing link to the periphery). 
Note that nodes 1 and 2 play a similar role and thus they receive the 
same score regardless of the navigation method at work. The differ- 
ence between the navigation schemes becomes more evident for the 
4 nodes in the periphery. RW assigns a similar rank to nodes 8 and 
6 and the couple 5 and 7, thus the two last positions of the RW rank 
are degenerated. However, in these two cases, these degeneracies are 
not supported in any similarity between the roles of the nodes. This 
problem is partially solved by PR which is able to unveil differences 
in the role of nodes 5 and 7 remaining, again, nodes 8 and 6 with the 
same score. Finally, QR solves the problem and assigns a different 
score to each of the 4 peripheral nodes. Importantly, QR considers 
node 5 as the most central element of the periphery, at variance with 
PR for which this position was occupied by node 7. 



PR since only nodes 6 and 8 remain with the same scores. 
However, the solution to the degeneracy of nodes 5 and 7 is 
not satisfactory enough as PR considers node 7 the most im- 
portant of the peripheral elements. A more reliable ranking 
is provided by the nonlocal nature of QR. First, the two de- 
generacies are finally distorted so that each of the 4 periph- 
eral nodes have different scores. Additionally, the hierarchy 
among the peripheral nodes is headed, as one would expect, 
by node 5 followed by node 7, whereas node 6 (a simple con- 
nector node between nodes 7 and 5) remains as the less im- 
portant node of the system. 
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E. QR in real networks 

To test the robustness of our previous findings it is con- 
venient to implement the QR navigation in real-world net- 
works in which more complex connectivity patterns appear 
In particular, we have considered a transportation network 
composed of the N = 500 major commercial airports in the 
United States and their connections via direct flights Isoll. In 
this case the key point is to detect those airports that play a 
central role in the complex organization of air transportation 
across the US. 

First we revisit the effects that the quantum-classical inter- 
play of QR causes on the time of convergence to a stationary 
ranking. In Fig. 3. a we have plotted the ratio t'^^/t^^ as 
a function of the classical-quantum interplay a. The curve 
confirms the speed up shown in Fig. 1 when a moderate quan- 
tum interplay is present in the QR navigation. In addition to 
this enhancement, we observed in Fig. 2 that the introduc- 
tion of quantumness produces relevant changes between the 
classifications obtained with PR and QR. For large networks, 
such as that of the US airports, the comparison between QR 
and PR is not as straightforward as for the simple graph of 
Fig. 2. However, a coarse-grained evaluation of the differ- 
ence between QR and PR can be obtained by means of the 
Kendall coefficient jsill . K, which takes a value K ^ 1 for 
perfect concordance between two rankings while K ~ for 
null agreement. In Fig. 3.b we have plotted the value of the 
Kendall coefficient between the classifications of QR and PR 
as a function of the quantum-classical interplay of QR, K{a). 
We observe that the values of K remain large for the values of 
a in which QR outperforms PR, but small differences appear 
as soon as quantum effects enter into play. 

Having quantified the differences that QR introduces in the 
final ranking, we now explore in more detail the nature of the 
changes. First we tackle the problem of degeneracies to test if, 
as shown in Fig. 2, QR provides a classification without less 
degenerate positions than PR. In Fig. 3.c we show the number 
of nodes occupying each of the positions contained in both 
PR and QR. This comparison reveals that, again, QR splits 
some of the degeneracies found in the classification generated 
by PR. This results is more obvious from the inset of the fig- 
ure as the classification of QR spans around 10 positions more 
than that obtained by PR. Thus, also in large networks, quan- 
tumness produces a finer discrimination of nodes centralities. 

Apart from solving degeneracies, the differences between 
QR and PR unveiled by the Kendall coefficient of both rank- 
ings raise the interesting question of which nodes are affected 
when introducing quantum effects into the PR navigation. In 
Figure 4. a we show the difference between the position of a 
node in the QR and the one assigned by PR as a function of 
the position in the first classification. The nodes ranked in the 
top part are not affected while most of the rises and falls are 
concentrated in the medium and lowest part of the ranking. 
This is, again, in agreement with the results found in the small 
graph in Figure 2. However, the case of the air transportation 
network is far more complex than that of the simple toy graph 
and we can gain more insight about the positive and negative 
effects that QR causes in the classification of lowly connected 
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FIG. 3. Quantitative differences between tlie classifications of QR 

and PR. In this figure we show the quantitative characterization of 
the perfomance of QR in a real network: the US air transportation 
system. In panel (a) we show the evolution of the ratio between the 
convergence times of QR and PR, t^^/t^^, as a function of the 
quantum-classical interplay, a. From this panel we observe, as in 
Fig. 1, that QR improves the convergence time of PR reaching the 
maximum outperformance when q°^* ~ 0.75. In panel (b) we com- 
pare the rankings obtained by PR and QR by means of their Kendall 
coefficient, K, which takes a value K = 1 for a perfect concordance 
between two classifications whereas K = for null agreement. As 
the panel shows, K decreases monotonously with a pointing out that 
the classification changes as quantumness is introduced in QR. Fi- 
nally, panel (c) represents the degeneracy found in each of the posi- 
tions of both QR and PR classifications. The fact that QR yields a 
ranking with more positions (see inset) implies that its rank is less 
degenerated than that of PR. 



nodes. 

In Figures 4.b-d we plot the US air transportation network 
by coloring each node according to the sign of the differ- 
ence between its QR and PR rankings. Those nodes in red 
have fallen in the QR (with respect to their position in the 
PR) whereas, for those nodes in green, QR has improved their 
PR score. Finally, those nodes being roughly at the same po- 
sition for both PR and QR are colored in yellow. First we 
notice that the structure of the network is governed by a set 
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FIG. 4. Qualitative differences between the classifications of QR and PR. Panel (a) shows, for each node in the US air transportation 
network, the difference between the position assigned by QR and that obtained via PR as a function of the first one. It is clear that for the 
most central nodes there are few differences between the position assigned by both navigation schemes. However, for nodes with moderate 
and low centrality the differences between QR and PR show up. In panel (b) we plot the US transportation network. The size of the nodes is 
proportional to their QR score whereas they are colored in green, yellow and red when the difference between QR and PR positions is positive, 
zero or negative respectively. In panel (c) and (d) we plot again the network but filtering those nodes with positive and negative values for the 
former difference respectively. To gain more insight, in panel (e) we plot the QR of each node as a function of the average QR of its nearest 
neighbors, QR"", divided by their average degree, fc"". The color of each point denotes the difference in the positions assigned by QR and 
PR. From this panel we observe that those points in green spans along the horizontal axis denoting that the nodes that are favored by QR are 
those connected to nodes with large QR but moderate connectivity [as it can be visualized from panels (b), (c) and (d)]. However, those points 
in red are localized in the bottom-left corner of the figure pointing out that those nodes punished by QR are connected to extremely large 
degree nodes (hubs) as, again, can be easily visualized in panels (b), (c) and (d). Finally, those nodes that are unaffected by QR (yellow points) 
correspond to large QR scores. 



(central part of the plot) of highly connected nodes (the hubs) 
whose degree is remarkably large. This is the set of large- 
score nodes whose elements do suffer any relevant change in 
their position. Additionally, we have a set of lowly connected 
nodes whose connections are mainly directed to the largest 
hubs. Most of these nodes are colored in red since they have 
been relegated to lower rank positions by QR. The fall of these 
nodes is balanced by the raise in QR of another set of lowly 
connected nodes (colored in green) located in the periphery 
of the network. In Figure 4.e. we analyze this effect in more 
detail. The nodes with the largest raise are those connected 
to neighbors with large QR and moderate degree, while the 
nodes that drop in the ranking are those that are connected to 
hubs distributing their influence amongst a large number of 
connections. Thus, QR ranks low degree nodes as a product 
of their global importance rather than to their local proximity 
to the most important elements of the system. 



III. DISCUSSION 

In summary, we have proposed a new way to navigate com- 
plex networks via quantum walkers. Our idea extends classi- 
cal navigation processes to the quantum domain, incorporat- 
ing coherences in the navigation through an otherwise classi- 
cal set of nodes. We have first introduced a framework for the 
walkers motion that combines the quantum unitary evolution 
with the iiTeversible character of the classical walk. This hy- 
brid navigation provides a unique stationary distribution for 
the occupation probabilities that allow to compute a reliable 
quantum rank. 

We have shown that the interplay between quantum and 
classical features in the random walk dynamics produces a 
significant speed up in the convergence to the stationary state 
of the walkers motion, thus decreasing the time needed to 
compute the final rank. We emphasize that QR can be read- 
ily implemented without the need of any quantum hardware, 
which does not exist yet, going beyond the current paradigm 
of quantum algorithms. 
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The introduction of quantum rules in the walkers also man- 
ifests the effect of quantum coherences in the properties of the 
rank itself. First we have shown that QR resolves undesired 
degeneracies found in ranks obtained with classical mean- 
ings. These degeneracies mainly affect those lowly connected 
nodes and it is within this subset where QR reveals more dif- 
ferences with its classical counterparts. Whereas those ele- 
ments occupying the top part of the classical rank remain as 
the most central ones for QR, the rest of elements are sub- 
jected to important variations. In particular, we have observed 
that those elements which are connected to moderately con- 
nected nodes increase their popularity with QR while those 
nodes connected to the hubs of the network decrease their 
importance. Thus, quantumness seem to punish the fact of 
being connected to influential nodes when this link becomes 
diluted by the large popularity of the hub. This punishment is 
balanced with the increase of the status of those nodes being 
connected to moderately central ones. Concisely, this result 
unveils the quantum solution to a recommendation dilemma: 
quantumness sets nodes recommended by locally important 
elements apart from those recommended by hubs which, in its 
turn, endorse a huge amount of leaves. 

As mentioned above, our algorithm does not rely on a quan- 
tum implementation. Therefore, it is complementary to recent 
studies on quantum complex networks focusing on the devel- 
opment of a quantum intemet ^3% . The backbone of these 
quantum networks is composed of photon emitters and re- 
ceivers with the aim of developing a new technology paradigm 
based on secure quantum cryptography protocols 1331. Obvi- 
ously, as the number of nodes and connections increases, the 
quantum information community meets complex networked 
structures. This synergy has lead to the emergence of an ex- 
citing new field in which classical complex phenomena such 
as percolation isil [35ll or the small-world effect lf36[ are 
revisited and addressed together with the new quantum ma- 
chinery underlying these networks. Here we have opened 
the door between quantum and complexity sciences from the 
other side. We have considered the actual network configu- 
rations and have applied a quantum perspective to the solu- 
tion of a classical complex network problem, trying to under- 
stand whether quantum mechanics may improve the solution 
to open problems of inherent complex nature Ii38i1 . 



Appendix A: Markovian Quantum Master equation and 
Random WalliS in Complex networlts 

Here, we discuss the theoretical framework for our algo- 
rithm. We introduce the markovian quantum master equa- 
tion (MQME) that generalizes its classical countenpart into 
the quantum domain. We demonstrate that for Google type 
transition matrices, the MQME is relaxing, i.e. it has always 
a unique stationary state. Finally we tackle the convergence 
times for reaching such stationary solution. 

To start with, we rewrite the master equation Q, for a clas- 
sical markovian random process: 



dr 



ijPj 



(Al) 



with 1\I = G — I, being G the Google matrix. The usual 
way of quantizing such random walk is by introducing the 
Shrodinger equation dt \ip) = —i/hH\ip), where the occupa- 
tion probabilities pi are replaced by the components of 
the wave function, such that. 



-iJ2^imi)m^ (A2) 



To finish the identification with iAli one sets {i\H\j) = Mij. 
At this point is mandatory to recall that H is an Hermitian 
operator, namely H = H\ implying that {i\H\i) = {j\H\i) 
(the matrix Al is real). Therefore the above procedure is re- 
stricted to undirected graphs, Mij = Mji. Similar problems 
arise with the discrete (coined) version of random walks ifisll . 
To solve the problem of directionality one can resort to the 
Szegedy recipe ll39ll . However, in any of the above cases, the 
evolution is unitary. In other words, classical-quantum iden- 
tification is not a straightforward task, classical Eg. (|Alb is 
irreversible while those quantum ones are reversible 12111 . 

There is a way of introducing both irreversibility and quan- 
tum coherences in a close way as it was recently pointed ll20ll . 
The idea is to work with Markovian master equations for 
the density matrix g rather than with the aforementioned 
Schrodinger equation. Being definite, the general theoretical 
framework is based in a time local master equation. 



dg 
'dt 



(A3) 
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Cg^-i{l-a)[H,g] (A4) 



where 7(i.j) > and are linear operators, acting on 

a finite dimensional linear (Banach) space: B{'H). The two- 
index notation (i, will be clear now. Equations (|A3l l and 
(IA4l i yield an irreversible dynamics where the directionality of 
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links can be readily implemented. In order to match classical 
and quantum random walks, we write the set of equations for 
the diagonal elements of the density matrix, g. Let us consider 
a = 1 in ( |A4] | and L(, = Then, 



(A5) 



In this case the equations for the non diagonal elements, gij 
{i 7^ j) become decoupled from the diagonal ones. By sim- 
ple inspection of dAll i and recalling qu = pi we identify 
7ij = Gij and both the classical equations (|2]i and the latter 
(|A5l l match. It is worth mentioning that, for a ~ 1, the non- 
diagonal elements decay exponentially in time, thus they do 
not contribute to the stationary state. Denoting g* the station- 
ary state, such that dt g* = 0, it is not hard to show (see below 
for a deeper discussion) that any initial state will asymptoti- 
cally approach to g*. Therefore in the discussed limit so far 
(a = 1) the stationary state has the non diagonal elements 
zero, g*j = i ^ j, while the diagonal elements g*^ = p* 
equal the classical random walk. 

The latter case is a trivial extension of the classical navi- 
gation process, without any quantum ingredient for the rank- 
ing. In general, we will deal with a 1. We choose for the 
Hamiltonian, Hij = 1 if the nodes i and j are connected and 
zero otherwise, in complete analogy to the identification used 
in quantum walks. It should be noted here that the hamilto- 
nian mixes both diagonal and non-diagonal elements making 
the dynamics non trivial and different from the classical one. 
Indeed, going to the limiting case of a one finds 



dt 



Qi'. 



-i ^2H,j{gj, - g,j) 



(A6) 



where the diagonal and non diagonal elements are coupled. 

1. Existence and uniqueness of the stationary solution 

We now prove a key result for our ranking: the existence 
and uniqueness of the stationary solution in the quantum ran- 
dom walk, g*. In order to do so, we start with the following 
definition: 

A master equation (23 with C given in A44I ) is relaxing if 
there exists a unique (steady) state g* such that any initial 
state go converges to g* for sufficiently large times. 

Besides we recall a Theorem shown by Spohn ll40ll : 

The master equation ( IA3I > with C given in ( IA4I ) is relaxing 
if ( i) the set j) } is self-adjoint ( the adjoint of every i(i.j), 
namely L^(^i jy inside the set) and (ii) the only operators 
commuting with all the are proportional to the identity. 

By putting all together, we state: 

(Theorem) A quantum random walk with master equation 
f |A3| ) and C given in ( IA4I ), where L(^i j-^ = \i){j \ and^ij = Gij 
being G the Google matrix d?]) is relaxing, i.e., it has always a 
stationary state and it is unique. 

(Proof) We first show that the set is self-adjoint. Given 



its self adjoint is L[.^^.j = = L 



We notice that the Google matrix links the nodes all-to-all. 
Therefore Gij ^ 0\/i,j providing that the adjoint operator 
is also in the set. Finally, be A an operator on B{'H), with 
decomposition: A = 'Y^aki\k){l\. Make the commutator with 



L 



(id) 



[A,L,,]=Y.'"'i[\k){l\A^){j\ 



k,l 

= Y.{au\l){j\-a,i\i){l\) 

+ (fli, - ajj)\i){j\ + aij\j){j\ - a*j\i){i\ 

since this conmutator should be zero for every i(i j) the only 
solution possible is that the non diagonal elements of 
(i 7^ j) vanish and its diagonal ones are equal, an = ajj. 
Therefore, A is proportional to the identity. Consequently we 
met the two conditions in the Spohn theorem, proving our re- 
sult. 



2. Approach to the stationary solution 



Being C time independent, the formal solution of ( IA4b is 



git) 



(A7) 



given any initial state g{to). It turns out that C can always be 
decomposed into a direct sum of Jordan formsj^: S~^CS = 
/:(°)©£(i)®...©£(-^), where. 



/Afc 1 

A;, 1 



V 







\ 





1 

A./ 



(AS) 



and C^^^ = (1), i.e. a 1 x 1 matrix with entrance 1, corre- 
sponding to the stationary solution (C relaxing) . Written in 
this way, the exponential yields e^* = 5(1) e (e^i*e'^i*) ® 
... © (e^^*e^^*)S'"\ with Nfe nilpotent matrices, so that 
the convergence to the unique stationary solution is guaran- 
teed since, limt^oo = c^' = S(l © © ... © 0)S"\ see 
ll22ll . The convergence ratio depends on how fast the blocks 
j,AfctpNfct approach to zero, so that are bounded by the largest 
eigenvalue (different from zero) Ai. This Jordan-Block evo- 
lution resembles the classical case, where the power method 
is widely used. Similar procedure could be used here by split- 
ting the evolution in discrete time steps: c'-* = ^m=ic'~^ 
with AI5 = t. In addition one could study the convergence 
by means of other Appendix, e.g. the proposed quantum adi- 
abatic algorithms 141] . 

We choose to integrate numerically the differential form 
and look for the convergence to the stationary solution. We 
define the convergence time r'^^ such that for t > t'^^ then 
\\git)~g*\\ <e, with the norm ||(?|| = TOft^I^P"^^. 
To compare both PR and QR we take the ratio 



where is the convergence time, integrating the classi- tion when = 1 see ( IA5I ). As initial condition we choose 
cal random walk, or equivalently the quantum master equa- = S^j/N. 
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